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' Abstract 

• We continue the analysis of perturbations in vector inflation. The dominant theme of this paper 

is the long wavelength limit of perturbations in small fields inflation and the controversial issue of its 
linear stability. We explain the nature of longitudinal modes, describe how they evolve, and show 
that they are not as harmful as it could seem at the first glance. On the other hand, the gravitational 
waves instability in large fields models is shown explicitly. It strongly limits a potential applicability 
of the recently proposed <5iV-type approach to vector inflationary perturbations. Finally, we expose 
a problem of an extra (gravitational) degree of freedom which appears whenever the vector fields are 
non-minimally coupled to gravity. 
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d . 1 Introduction 

, Recently it became evident that higher spin fields can source inflationary expansion of the Universe 

;> ■ [U [2l [3]. At the level of background FRW dynamics, higher spin models go almost identically to the 

scalar infiaton case. But they also provide a possible account for a large scale anisotropy [T] which is 
currently of a great phenomcnological importance. It motivated many independent researchers to study 
cosmological perturbations in higher spin inflationary scenarios [U[5J[51[7]- However, the problem appears 
to be very complicated due to non-trivial couplings of different types of perturbations (scalar, vector, 
tensor) to each other even at the linear order [5], with a possible exception of a 3-form infiaton field 
■ [2j [5] being dua{3 to some peculiar scalar [2l [7] . 



| In this paper we focus on vector inflation with the action [TJ [U [5] 
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where /(„) = —A^A?-. and Fj$ = V^Au — VvA^; the vector fields are supposed to be randomly 
oriented so that the background metric is approximately isotropic with linear perturbations given by 

ds 2 = (1 + 20) dt 2 + 2a(t)V l dtdx l - a 2 (t) ((1 - 2ip) 8 lJ - %) dx l dx j , (2) 

V\ = 0, h\ = 0, = 0. The mass-term inflation corresponds to V = —^A^A^ — ^1. For every 
single infiaton field the energy-momentum tensor takes the form 



+ i (m - l -51R^ A 7 Ay + \ (S$D - V Q V ) A 7 A, (3) 

1 Actually, the duality is flawed by higher time derivatives in the dual action. These higher derivative effects seem to be 
intrinsically non-linear and do not spoil the linear analysis [7]- Probably, they can be consistently eliminated in perturbation 
theory [10| . But it is also quite possible that this problem reveals an actual extra degree of freedom which is discussed for 
vector fields in Section 5. 
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which is by itself quite capable of giving the general feeling for why the perturbation theory is so messy. 
A crucial simplification can be achieved assuming that inflation is driven by small vector fields, NB 2 <C f , 
B = . Anyway, this assumption is almost unavoidable if one wants to ensure stability of gravitational 
waves f4]. In Section 2 we accept it and consider the perturbations for very small values of inflaton fields. 

After that we proceed with elucidating some tricky aspects of vector inflation. Namely, in Section 
3 we give a thorough analysis of longitudinal modes and associated stability problems [HI [12] , and in 
Section 4 we show the gravitational instability of large fields vector inflation explicitly, for it has recently 
been doubted in [5j. After that, in Section 5 we report a new problem of vector inflation concerning the 
number of degrees of freedom. And in Section 6 we conclude. 



2 Linear perturbations in small fields inflation 

The full (and horrible) set of linear perturbation equations can be found in [5]. Fortunately, we need only 
a few rudiments of the general perturbation analysis. Every term in ([3]) can be varied easily, although the 
whole expression becomes very bulky. It is clear that scalar, vector and tensor perturbations mix with 
each other because we can contract a background vector Bi with a perturbation. For example, using the 
metric ((2} we have SB 2 = 2B l 5B l + 2B 2 ip + h^B^j. And in another linear relation, SA° = SA + Vj-Bi, 
the quantity ViBi is also a scalar. From A = aB and H = - it follows that the leading contributions 
of most of the terms in ([3]) are proportional to H 2 B 2 , with indices contracted or not. However, at the 
background level the largest terms cancel each other (with ~ -7= accuracy for N random fields and 
exactly for the fine-tuned case of a triad), and the background dynamics coincides with that of scalar 
N-flation pQ. As for the fluctuations, one has of course to take H 2 BSB terms into account, and we will 
see below that they are important. 

While analysing the possible perturbations in §3§ step by step, it is rather tempting to conclude that 
the curvature perturbations in vector inflation with a large number of fields are ridiculously small. For 
example, one could argue that any terms in Too of the form of H 2 ^2 AjSA^ are statistically suppressed as 

N — > 00 because the fluctuations 5 A can have arbitrary directions. However, it is not a reliable argument. 
Indeed, we definitely want to keep the Hubble constant intact in the course of the limiting procedure (or 
at least the Hubble rate should not diverge) . In the mass-term inflation it means that A oc 1/ v~N, and 
despite the 1/VN statistical suppression the length fluctuation term H 2 ^ AjSA^ has no scaling with N. 
It shows that one has to look for some other approximations. 

Note that this type of naive argument, if it were only correct, would also suppress the perturbations 
in scalar N-flation in contradiction to the general statement of [13] . And the reason for which it does not 
actually happen is precisely the same as for vectors. Namely, for V = ^- Y^, 4% we S et V x N^cj) 2 and 
SV cx \/Nm 2 (j)(S(f>) where (S(j)) is a typical magnitude (variance) of fluctuations. And then for the relative 
magnitude of perturbations we have ^ cx -7= ^^r^- There is the \ZjV factor in denominator. However, 
recall that H 2 cx m 2 N(f) 2 and hence <j> cx Let's parametrize the wavelengths A by K = (XH)^ 1 

so that k = 1 at the horizon scale. Then we have 5(j) K cx kH cx K\fNmcj) and finally cx nm. The 
magnitude of perturbations depends solely on the inflaton mass. 

In vector inflation all types of perturbations are mixed, and it is particularly interesting to evaluate the 
effect of mixing with gravitational waves due to H 2 ^2 BiSBj contributions to the linear fluctuations of 
the stress tensor ([3]). If we assume the fluctuations SB are random (e.g. their directions are not correlated 
with the background direction of the field) then this term has the usual — j= suppression. So that the 

effect should be of order ^/NH 2 B(SB). For the mass-term inflation we get 5B K cx kH cx KmB-^N, 
and the relative weight of this perturbation >/~N {B5B) cx nmB 2 N is hugcl when we start inflation at 
B ~ N^ 1 / 4 . Some perturbations are completely out of control at the onset of large fields inflation. 
It is understandable because the whole story of the background dynamics has emerged from statistical 
cancellation of the leading terms cx H 2 B 2 in the energy-momentum tensor, and inflation is initiated when 
anisotropic corrections are of order one (and actually badly unstable, see Section 4). In new inflation 
the quantity \^NB5B cx kH^/NB can be made relatively small. This term serves as a source for gravity 
waves provided by inflatons fluctuations. One can control its magnitude by varying the form of potential 



2 This can raise some doubts about the instability of gravitational waves in such inflation as it was deduced neglecting 
these terms 0. Intuitively it's hard to believe that some external force can neutralize the effect of the large tachyonic mass. 
And in fact, this intuition works well, see Section 4. 
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and the number of fields in order to produce a desired amount of tensor perturbations which is usualy a 
problem in new inflation (see also [14] ). 

It is very natural that, under suitable conditions in the small fields limit, the gravity waves disentangle 
from the other modes. After all, the linear mixing of different perturbations occurs due to the presence of 
(a random set of) preferred directions. But the vectors which represent the preferred directions become 
very small and not very important, therefore the mixing of modes is notably weak in this limit [5]. Hence 
we can approach the dynamics by considering background values of the inflatons as small quantities. 
(Although they should be larger than the perturbations, of course.) For the equations of motion of the 
vector fields it means that in the first approximation we are to consider the fixed background geometry, 
i.e. neglect the gravitational backreaction. Indeed, all other terms in the first-order equations of motion 
contain variations of metric multiplied by background vector fields, see also [5]. In this limit we substitute 
the non-perturbed FRW metric (with a fixed but otherwise arbitrary time-dependence of the scale factor) 
in the action JI} and get for the vector fields pQ: 

- ^AA a + \2Vj + |^)Ao + ~ (diBi + HdiB^j = 0, (4) 

Bi + 3HB t + 2VJB, - -^ABi - - (ftio + Hd t A ^ + ^ (djBj) = 0. (5) 
And there is also a consistency condition 

V|1 hvj + | \ A» = 0. 

Note that for H = const one gets a de Sitter background which amounts to setting the slow-roll parameters 
to zero. 

From now on we would be interested in only the leading behaviour of the vector fields and therefore 
assume the geometry to be a pure de Sitter (so that the scalar curvature is constant in time) and neglect 
the variation of effective mass (it would contain Vn and extra powers of A). Then we have V p-A? = 0, 
or more explicitly 

A + 3HA - -d^ = 0. (6) 
a 

It allows to simplify the equation 

1 2H 

Bi + 3HBi + 2V I B l - -tAB, + = 0, (7) 

a z a 

and 2Vj can be substituted by a constant mass. We see that in the long wavelength limit all component^] 
of Bi evolve exactly like scalar fileds with the mass m 2 — 2Vj. The dominant mode slowly rolls, while 
the other one fastly decays. And the condition © (or, even better, the constraint equation ^j) shows 
that Aq decays exponentially in physical time as it was correctly stated in [5]. However, it is not true 
for the spatial longitudinal mode itself, as it behaves identically to the transverse ones when A — > oo. It 
is just its contribution to the constraint equation © what goes to zero. And actually, it goes to zero as 
- (this dependence represents the stretching of waves) so that A cx - also (more precisely, a little bit 
different from that due to the slow motion of the field B). Moreover, the mode which could be erroneously 
deduced from © naively setting the spatial derivatives to zero, namely Aq cx ^j, is absolutely fake. The 
consistency condition is necessary but not sufficient for a vector field to satisfy the whole system of 
equations of motion. From ([4]) we see that if there is no longitudinal mode, then Aq = exactly. The 
superhorizon analysis is always notably subtle in that it is not reliable to neglect the spatial derivatives 
for evaluation of subleading quantities. 

For the Einstein equations we can use the same philosophy and drop the terms which contain two 
powers of B : like B 2 ip; then, for example, S(B 2 ) would be given just by 2BSB. If we could also neglect 
the gravitational waves and consider the anisotropy as only a small correction, then the number of e-folds 
in any patch of the Universe would be well-defined and depend only on variations of Bi (recall that A 
decays). In this limit and under these assumptions one can use the (W-approach of [6J. However, if the 
anisotropics grow as fast as for the mass-term inflation then it would be a tricky business even to speak 
about the number of e-folds. In the limit of small fields the accuracy of <5Af-approach is determined by the 
strength of mixing with tensor modes. And as the latter can not be much greater than the tensor-to-scalar 

3 See the next Section for peculiar properties of the longitudinal component which arc missed in this limit. 
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ratio (for not to produce too much of gravity waves), the SN formalism is applicable to at least a few 
percent accuracy for viable models, and we can be sure in the nearly flat primodial spectrum of curvature 
perturbations if the gravitational waves are tame. Moreover, T 00 (see [5]) depends then only on the length 
fluctuations BjdBj (and their time derivatives) which explains why the results of [6] are identical to that 
for scalar N-flation. All the relevant variations in the slow-roll regime come from changing the scalar 
argument A 2 of the potential. In general, precision of these results can be roughly estimated as above, 
and beyond this accuracy the perturbations differ from that of scalar inflation, although the detailed 
predictions are model dependent and presumably require numerical methods. 

Finally, we should mention the vector perturbations of the metric. They usually decay as fast as 
(see for example [15] ) obeying the conservation law of angular momentum. The vector inflation is no 
exception, although the decay could be a bit slower due to vector perturbations of the energy-momentum 
tensor. In absence of vectorial stress tensor perturbations the spatial part of Einstein equations would 
give the usual decay. But the right hand side of it contains in vector inflation a source for V. The source 
represents the vortical excitations of vector fields which carry non-zero angular momentum. Therefore, 
the vector inflatons can produce some vorticity. But apparently, they can not make any kind of vortical 
instability at least in small fields models because, neglecting the gravitational backreaction, the inflatons 
evolve in FRW space-time and by themselves obey the angular momentum conservation, so that the 
source for V also decays during expansion in much the same way as a rotating platform slows down if 
some mass moves on it in the radial direction. The backreaction would allow the vector fields to gain some 
vorticity by transfering the opposite angular momentum to the gravitational field as a recoil. Therefore 
we expect that for vector inflation the vector perturbations decay can be a bit slower than the usual 
^j-law, but not much slower for small fields models. It is hard to make more precise statements on V 
due to the aforementioned problems of working with subleading quantities in the long wavelength limit. 



3 The problem of longitudinal components 

At a finite wavelength, the equation (fj| and the scalar part of equation ([7]) describe the evolution of 
temporal and longitudinal components of vector fields. Note that (T5)) is a constraint equation, and 
therefore we consider only one degree of freedom per every vector field (see, however, Section 5) which is 
actually suspected to be badly unstable for tachyonic masses [TTJ [12] . It was claimed to be a ghost which 
would make any attempt to pursue a quantum field theory nonsensical, at least in the usual perturbative 
formulation. In this Section we discuss how to live with it. 



3.1 Vector fields in flat space-time 

To set up the stage, we first describe a tachyonic (m 2 < 0) vector field in Minkowski space. We set 
R = 0, H = 0, B = A and a = 1 in the equation of motion Q and get (□ + m 2 ) A = for the 
spatial components of the vector field. It is just the tachyonic Klein-Gordon equation with no apparent 
problems. The problem comes from equation (— A +m 2 ) Aq + diAi = 0. Fourier decomposition in 
spatial modes cx e lkiXi shows that the temporal component 

A ° = -12X^2 8 

K + m l 

diverges at fc 2 = — m 2 > 0. At this wavelength the longitudinal mode is not permitted (but Aq is 
arbitrary) and in the neighbourhood of fc 2 = — m 2 it involves very large values of Aq. As long as the 
tachyonic vector field is free, one can quantize just three independent fields Ai and forget about unphysical 
variable Aq but this approach is not suitable for turning on Lorentz invariant interactions. 

In principle, the above result looks like a kind of resonant behaviour which one could try to kill by 
a suitable non-linearity. For example, it is tempting to construct a potential which would give a normal 
mass at large values of Aq. Unfortunately, given the Lorentz invariance, the best we can do is to make it 
so for large positive A^A^. Then the differential operator — A +2Vj would still have zero modes at the 
length scale of |ro e //| 1 for negative values of A 2 where the field is tachyonic. The temporal component, 
given by ([8]) , would have to be very large around these modes which is compatible with negative A 2 and 

is even largeiQ. Thus, the best we can do is to shift the problem 



tachyonic mass if the spatial length 



.4 



4 Of course, Ao is by itself proportional to the length of longitudinal mode but (in dimensions greater than 1 + 1) we can 
always play with the lengths of transverse components to make A 2 negative. 
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to the region of very large fields, ft is not possible to completely exorcise the instability with modified 
potential, at least without abandoning Lorentz invariance or introducing new degrees of freedom. The 
Hamiltonian analysis also confirms this result. 

3.2 Hamiltonian analysis 

As usual, the conjugate momenta are defined as 7r M = -^j- and given by 7Tj = F oi and n = (primary 
constraint). For the mass-term potential the Hamiltonian density reads 

I 1 771? 

H = -tt? + (diAo)^ + -F^ + — (A 2 - A 2 ) . 

Commutation of the primary constraint tt with the Hamiltonian gives the secondary constraint^ d^i + 
m 2 Ao = which allows for elimination of unphysical variable A$: 

H = y + i {d ^ )2 + \ FijFij + 

If m 2 > then all the terms are positive, while for tachyonic masses the second and the fourth terms 
are negative. The latter is just the ordinary tachyonic potential but the former contains momenta and 
therefore looks like a ghost. However, the secondary constraint shows that for the physical states it equals 
just to rj ^Al: 




There is no problem in the infrared because A§ — > if k — > 0, and no problem in the ultraviolet because 
Aq cx 4^ cx 0(1) when k — > oo so that another term jFijFij cx k 2 A 2 wins the game. It is not a genuine 
ghost as it presents a problem only in a limited range of k 2 near — m 2 in a sharp contrast with real ghosts 
which go worse and worse in the ultraviolet. 

For a general potential the secondary constraint is 9j7Tj + 2VjAq = 0, and in terms of physical 
variables the Hamiltonian becomes quite complicated since we have to use Aq = —i^r~ m the argument 
of V. However, we can go the opposite way again and exclude diirt: 

n= 1 r 2 + 1 -F ijFij + ^- + v(i). 

Now one can easily convince himself that it would be rather difficult to make this Hamiltonian positive 
definite because at some point we need to pass from Vj < to Vj > 0. 

So, we do not understand the theory completely. But it is not a fatal problem since we do not 
necessarily have an infinite phase space volume in the ultraviolet for the large Aq instability. And the 
crucial property of inflationary space-times is that the wavelengths are changing with time and every 
mode dwells at a singular point k 2 = —m 2 only for a single instant of time. Now we proceed to show 
that it solves the instability problem at least at the level of classical equations of motion^ using the test 
field approximation in the simplest case of de Sitter geometry. 



3.3 Longitudinal mode in vector inflation 

For a longitudinal mode (kiBi = kB) we use the constraint (f4]) to find the closed form of the equation of 
motion in de Sitter space (f = -2H 2 ): 

.. / 2H4 \ ■ (k 2 2 2if 2 4 \ 

B+13H+ -T3 £ } B + + m 2 + -r- 2 2! S = 0. 9 

^ £+ m 2-2H 2 J \a 2 ML +m 2_2H 2 J 

At short wavelengths there is no big difference between transverse and longitudinal modes. The problem 
appears when the wavelength is near the value for which |r = 2H 2 — m 2 . It's time to pay for introducing 
the tachyonic effective mass ~2H 2 . We need to understand the properties of solutions in this region. 
Let's neglect the mass m of the inflaton for simplicity as it much smaller than H in the slow roll regime, 

5 As one can easily check, it coincides exactly with J3J. 

6 In [6] it is argued that even the quantum theory makes a good sense in this situation. Of course, these arguments can 
be taken seriously only after we have checked that we don't really deal with a ghost. 
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and pick up a wave which was well under the horizon with K = AH 2 at some instant of time to ■ Then 
^ = AH 2 e~ 2H ^~ to) and © takes the form: 

(2 - e 2H ^ to) ^j B + (lOH - ZHe 2H{t - to) ^j B + %H 2 e- 2H{t - to) B = 0. 

With a new time variable r = 2H(t — to) — hi2 we get 

2(l-e T )B + (5-3e T )B + e- T B = 0. (10) 

The critical point of crossing the singularity is at r = 0. The coefficient in front of B vanishes at this 
point and all trajectories are tangent there to a one-parameter family of curves B = — ^. This behaviour 

is stable because if B ^ — -j at small r < then the second derivative B ~ — —2- has an appropriate sign 
to correct the trajectory. 

Our task is to find a two-parametric family of solutions for equation (|10| . From the previous analysis 
it is clear that one possible solution contains no longitudial mode at all when ^2 = 2H 2 . We take 
B(t = 0) = B(t = 0) = and construct the solution in the form of power series B = t 2 + C n r n . 

n>3 

The first term solves the equation up to 0(T 2 )-terms, and O(r 3 )-corrections give C3 = — I and so on. 
This solution is rather smooth around the problematic point due to absence of longitudinal mode, it 
corresponds to Aq = at t = 0. 

The second solution is more interesting. We write it down as £> = 1 — ^ + D n T n which gives 

n>2 

finite values of Aq due to cancellation of two first-order zeros in (JSj) . The first two terms solve ([TO]) at 
the level of 0(t), then D2 is undetermined because at the level of 0(t 2 ) it solves the equation by itself, 
thus one can take D2 = and proceed with D%. The full two-parameter family of solutions is given by 
B = a — + /3t 2 + C n T n where a and (3 are arbitrary constants while C n -s should be determined 

n>3 

one by one in terms of a and f3. 

In order to understand the properties of the second solution around r = we note that, unlike for 
the first one, the main players for it are B and B. If we neglect the B-term in (fT0|) . then the resulting 
equation can be solved explicitly in elementary functions. We present the solution in terms of the physical 
time: 

One can easily check that at the moment of time t — to = Inl we have B = —BH as required. It is also 
straightforward to estimate the value of Aq at the same time: 

(B + HB)e H< - t - t ^ C 
~ (2-e 2ff ( t -*o))i7 ~ 20 

which is clearly not too large (although it is the value of a ^ fraction), and no catastrophe happens. The 
t — > 00 asymptotic is completely stable too, the amplitude goes to a constant value (because the inflaton 
mass was neglected) of the same order of magnitude as it was when we picked it up at t — to- (Of course, 
one should properly rearrange the signs for t > to + In ^.) The solutions can be evolved through the 
dangerous point smoothly. 

One could make the above analysis without neglecting the mass in ([7]). It would not change the 
qualitative results, but the coefficients would not be so nice. The passage through the singular point 
(with zero coefficient in front of B) gets shifted in time, but the signs of all the terms around this 
point are the same as before. Instead of the B = —HB behaviour at the crossing point, we would find 
B = —H ^1 + jj^j B, and in the infinite time limit the longitudinal mode would exhibit the slow roll 

evolution B ~ e~^H l . Of course, when taking into account the inflaton mass, one should better use the 
actual FRW-metric instead of the pure de Sitter one, see [IB] . 

There is no visible signature for the linear theory breakdown. Note however that the amplitude 
changes from its initial value to zero, and then to values of opposite sign at a time scale of order of one 
c-fold. Of course, this is just the amplitude of some Fourier modes around the horizon, so that it does 
not imply anything catastrophic, but in general the time variations of A 2 in the argument of potential 
term can be somehow more important than usually. It means that non-linearity of potential (changing of 
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effective mass) could play some role, and therefore considerable non-Gaussianities in primodial spectrum 
could be produced. 

When the present article was in preparation, another work has appeared on arXiv, namely [16j . in 
which the same problem is analysed. The Authors of [16] also came to the conclusion that the linear 
longitudinal modes safely pass through the point of k 2 + M 2 ^ — 0, at least in the case of a test field in 
FRW- Universe which is dubbed the zero vector vev in their work. They also argue that in other situations 
it is not the case: they claim that even a single vector field with non- vanishing vev in Bianchi I space-time 
develops a singularity of linear modes at this point. It is quite contrary to our general intuition which 
assumes that inclusion of the metric perturbations in (fT0|) should not make the transition harder, recall 
that it is basically governed by the sign of the coefficient in front of B. We would even expect that the 
transition should be softer in presence of metric perturbations since not only different modes but also 
different parts of the space would pass through this point at different instants of time. However, the linear 
instability was inferred in [16j from numerical simulations. The Authors do not expain their numerics 
explicitly, but one can safely guess that they evolve the longitudinal mode from deep inside the horizon 
directly solving the equations of motion by standard methods which involve explicit determination of the 
second time derivatives of the physical variables via the other terms in equations. It basically amounts to 
solving our equation (fTU)) for B and using this value at each step to construct the numerical solution. We 
suspect that, probably, this numerics fails to give a meaningful answer for a very complicated system of 
equations near the point at which the expression for the highest derivative terms is singular. On the other 
hand, one should also take care about a possible anisotropic instability in the system (see the Section 4), 
which indeed could change a lot. 

Note also that one more problem with linear analysis was pointed out in |16j . The linear solution 
diverges when M^f* = m 2 + = 0, approximately at the end of inflation. We want to explain here the 

reason for that, which is actually quite simpleQ. The consistency equation y M ^M^^A^j = reduces at 

this point to dM ^ t !1 A a = 0. And it clearly can not be satisfied with general initial conditions. This is a real 
problem which occurs only once at the exit from inflation. And again, it signals the lack of fundamental 
understanding of the nature of the vector fields in the model. At this (and only this) instant of time the 
gauge freedom exists which makes the number of degrees of freedom ill-defined. The longitudinal mode 
becomes infinite exactly at the point at which it fails to be physical. 

4 Instability of large fields inflation 

In Ref. [4J it was shown that large fields vector inflation is badly unstable with respect to gravitational 
waves, i.e. small anisotropics grow in it much too fast. It was argued in [6] that this conclusion is just an 
artifact of the linear approximation technique in the Jordan frame, while in the Einstein frame everything 
should be stable because the interaction of gravity with matter has the normal form there. However, we 
would like to remind that, roughly speaking, only two thirds of the instability reported in [4] came from 
the |r A 2 term in the action, while a remaining one third came from the kinetic term of the vector fields 
which is conformally invariant. Of course, for this instability to develop we had to stabilize the fields B 
(to ensure the slow roll), but once it is done, the instability is there in both frames. In this Section we 
analyse the behaviour of anisotropies in a manner which is somewhat closer to <5A-formalism. 

We assume that in a separate patch of the Universe the Hubble rate experienced a sudden jump in 
only one direction, H z = H + h. The metric would be the axially symmetric Bianchi I type 

ds 2 = dt 2 - a 2 (t){dx 2 +dy 2 )- b 2 {t)dz 2 

with the new components of the Einstein tensor: G% = H 2 + 2H a H b = 3H 2 + 2Hh, G% = = 
H a + H b + H 2 + H 2 + H a H b = 2H + 3H 2 + h + ZHh + K 2 and G z z = 2H + 3H 2 . If we also assume 
that the energy-momentum tensor remains isotropic, then subtracting G z z from G x x gives the equation 
h + 3Hh + h 2 = which shows that anisotropies are being washed out (as they should be according to 
the general theorem of Wald, [U]). If H w const then h oc for small h. One could worry about large 
negative jumps which could become stable at h —3H, but it would already correspond to a phantom 
matter with energy density —3H 2 . And in any case, large anisotropies are not generally expected to 
decay since, for example, Kasner solutions are known. 

We claim that in vector inflation non-zero h would render the energy-momentum tensor anisotropic 
too. In the homogeneous limit and with the slow-roll assumption, the main contribution to anisotropic 

7 See also the last Hamiltonian of the Section 3.2. 
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part of T? comes from fAiAi, -F 0i F°i and \G\A 2 terms. The latter one rcnormalizes the gravitational 
constant, and we forget about it for a moment. But the first two terms give important contributions to 
T£ — T* . With a natural definition B z = ^ff- we get for them 

-Ef(^-^ 2 )-™ 2 ((^>-(^)), 

J2(H 2 X B 2 X - H 2 B 2 ) « NH 2 ({Bl) - (B 2 Z )) - 2NHh^- 

if both h and (B 2 ) — (B 2 ) are small. And we also have to take the subleading contribution to —Fq^F ^ 
into account, namely 2^2(H X B X B X — H Z B Z B Z ). It is easy to derive the equations of motion of vector 
fields: 

B x + (2H a + H b )B x + (m 2 + ^+H a + Hi + H a H^j B x = 0, 

B z + (2H a + H b )B z +(m 2 + ?-+H b + 2H a H b ^j B z = 0, 
which show that in the slow roll regime we have 

a _a „ H b -H a + H a H b -H 2 h + Hh 
2H a + H b ~ 3H 

for a pair of identical fields in directions of x and z axes. And therefore 

2 J2(H X B X B X - H Z B Z B Z ) w 2HBj ■ ^§^B w ^-(h + Hh). 

Note that at the moment of turning on the fluctuation of the Hubble rate we do not change the value of 
the scale factor itself, therefore we may keep both A and B independent of the spatial direction of the 
field. Then the anisotropy of T/ initially equals to -2NHh^- + 2N 9 B (h + Hh). Now we have to recall 
that the \G\A 2 term effectively divides the gravitational constant by 1 + and the resulting effective 
anisotropy would be —4Hh + |(A + Hh) in the large fields limit. It has a drastic effect of making the 
linear equation for h unstable: —\h + ^-Hh = 0. 

At the same time the anisotropy of matter distribution grows according to vector fields equations of 
motion: 

| ■ 3NH> « S *> - m ^.2B. h -±^B - (A + EH) , 

and in this case tends to stabilize the fluctuation. However, if we differentiate the anisotropic part of 
Einstein equations with respect to time: 



d I 1 l 17 TTiA 1 NB \ 2HNB 2 (■ x 

we get a simple approximate equation h — 5Hh + \2H 2 h = with initial condition h(0) = 17Hh(Q). Its 
solution exhibits a fast exponential growth with oscillations. 

Admittedly, the geometric effects which overturned the sign of h can be changed by, for example, 
taking anisotropic distributions of the vector fields. But even if one can force h to decay initially, 
then after a small period of time (before the perturbation could finally decay) one generically gets a 
considerable amount of anisotropy in the term which does not depend on the current value of h, namely 
3NH 2 ((B 2 ) — (B^ )) . This conclusion is robust. It can either prevent the Hubble rate jump from stopping 
at zero value and make it growing with the opposite sign, or change the sign of h before the fluctuation 
could reach the zero, depending on the coefficients in other terms. In any case, a small anisotropic 
fluctuation in the Hubble law is badly unstable for the large fields vector inflation. 



5 An extra degree of freedom 

Now we want to report on a new problem with vector inflation. A naive counting of independent propa- 
gating degrees of freedom would give the number of 3./V+2 for vector fields and one graviton. Equations 
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also contain the second time derivatives of tp. But it is usually not a dynamical quantity since the tempo- 
ral component of the Einstein equations contains no second time derivative^] and can be (and should be) 
regarded as a constraint which determines ip completely. (It actually corresponds to Newtonian potential 
defined by distribution of masses in non-relativistic limit.) However, for vector inflation it is not the case 
beyond the linear approximation since T 00 has a term with RAq and the scalar curvature depends on ip. 
On the other hand, the field equation (fj)) is also no longer a constraint due to precisely the same reason. 
And it is not clear how to find those linear combinations of equations which would give constraints in 
the Jordan frame. 

To find out the actual number of degrees of freedom we need to perform a conformal transformation 
to the Einstein frame: g^ v = e 2p g p „, R — e~ 2p (R — 6Dp — 6(Vp) 2 ), = A^' with the conformal 

weight 

P=^(l + it '(»))■ 

It transforms the original action tc^| 

s= /<W^( - f-v /Ir -v" r --v- VtC'C/.'/'V''! i" 1 



where 

and the new potential is given by 



\ m / n \ \ p 



(p) J • 7 («) 

so that the vector fields are intricately interacting now. This potential looks scary, and one can think of 

2 

the special case V = for which it is just 



Anyway, our discussion does not depend on a potential. It refers to any non-minimally coupled vector 
fields (and should be suitable for higher forms too). 

We first consider the new action (fTT|) with N = 1 in a flat (Minkowski) space-time. Canonical 

momenta are given by ttq — AAgip^ (^AqAq — AjA^j and TTi = —AAiip 2 j (^AqAq — AjAjj + A\ — OiAq. 

This system of linear differential equations for velocities has a determinant 

-{A^fAl (AlAl + AlAj + AjAl) 

which is not zero if Aq ^ and therefore is solvable under this assumption, and no constraints are there. 
Thus, the vector field acquires a fourth degree of freedom whenever the longitudinal mode is on. The 
reason is easy to understand analysing the equation of motion 

dip dV 
2A v -^(n<p)+V fi F» v + 2A»— = 0. (12) 

Its temporal component is not a constraint due to the time derivatives in □<£>. 

For an arbitrary number of fields, evaluation of the determinant is a tedious procedure. 3iV + 1 
degrees of freedom are guaranteed, but a priori there could be even more of them, up to AN. However, 



8 One notable exception is a treatment by Weinberg |18| who casts the equations into a form = T M „ — j Sm^^o* which 
contains time derivatives in the temporal component too. (Weinberg has also an overall minus sign in the right hand side 
due to the opposite sign convention for the Ricci tensor.) Of course, a proper linear combination of these equations restores 
the constraint. 

9 In the slow roll regime the conformal factor changes slowly with time, so that the scalar quantities A 2 and F 2 remain 
almost frozen in both frames. The new kinetic term (VV 3 ) 2 pl avs the role of non-minimal coupling ensuring the slow roll 
conditions. 
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we can check that there is only one extra degree of freedom at any value of N by uncovering N — 1 
constraints. It is easy to do. Suppose, we have two fields with equations of motion (fT5|) . Let's multiply 
the temporal component of the first equation by 2^4" 2 ^ , the temporal component of the second 

equation by 2A9js gj^- and subtract them from each other. Time derivatives are cancelled, and we get 
a constraint. Then we can play the same game with other fields. Clearly, there are N — 1 independent 
relations of this sort. 

We see that whenever at least one of the fields has Aq ^ 0, an extra degree of freedom turns on 
and completely decouples again in absence of longitudinal modes. It can't be seen at all at the level of 
background dynamics. Although this mode should not play any role in the linear perturbation analysis 
around the state with Aq = 0, it means that at the fundamental level the theory is not very well defined. 
The problematic degree of freedom is essentially gravitational because it is one for all the fields and stems 
from non-minimal coupling term due to time derivatives in R. Perhaps, it could be regulated by explicitly 
turning it on even in the Aq — > limit, for example with a small i? 2 -correction to Einstein gravity. 

6 Conclusions 

An interesting possibility of driving inflation with higher spin fields is nowadays being emergent. It 
appears that this idea can be perfectly realised at the background level with almost no complications 
[TJ [51 |2j and with a benefit of a natural large scale anisotropy suggested by some recent observations [H] . 
But already at the level of linear perturbations the full set of equations of motion [5] becomes almost 
analytically untractable (with a possible exception of 3-form inflation [2j [7l 13 E] ) ■ And also some problems 
have been reported both before [HJ [12l [TBI [20] and now, in this paper, see Section 5. However, we have 
shown that although there is obviously a serious lack of understanding the nature of vector inflatons at 
the fundamental level (let alone the problem of UV-completion) , nevertheless one can consistently use 
small fields inflation as an effective theory of the early cosmological evolution. 

Very recently, a first calculation of non-Gaussianities for inflation with vector fields has appeared, see 
|21j . It is claimed that vector fields can produce a high level of non-Gaussianity. The analysis in [21] 
is performed in the usual JTV-formalism of [6] with all its potential shortcomings, but the possibility is 
nevertheless very interesting. Note also that among the things which could not been taken into account in 
|21j is the tricky evolution of longitudinal modes. We have seen above that this evolution is very peculiar, 
and therefore a general intuition would imply that it also can make a contribution to non-Gaussian 
features. Of course, this issue needs a further investigation. 

We conclude that there are many unresolved fundamental issues about vector (and higher spin) 
inflation, but at the level of effective description it is a viable candidate for the theory of inflationary epoch. 
Moreover, it could provide interesting insights into quite a few problems of current phenomenological 
interest. A great deal of further progress remains to be done and, from our perspective, worth to be done 
in the field of understanding vector inflation. 

The Author is grateful to Vitaly Vanchurin, Viatcheslav Mukhanov, David Lyth, and Cristiano Ger- 
mani for very useful discussions. This work was supported in part by the Cluster of Excellence EXC 153 
"Origin and Structure of the Universe" . 
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